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Overview: 


This  report  contains  three  parts.  First,  a  list  of  archival  journal  publications,  conference 
proceedings,  and  abstracts  for  oral  presentations  are  listed.  It  is  noted  that  the  research  reported 
in  these  papers  and  presentations  was  performed  in  conjunction  with  “Generation  of  Synthetic 
SAS  Data  for  Targets  near  the  Seafloor:  Propagation  Component,”  (ONR  Award  N00014-07-G- 
0557/0023).  Second,  the  fast  ray  model  for  target  scattering  in  a  homogeneous  waveguide  has 
yet  to  be  published  in  an  archival  journal,  but  an  overview  has  been  given  in  the  conference 
proceedings  for  the  First  International  Conference  and  Exhibition  on  Underwater  Acoustics, 
Corfu,  Greece.  As  this  proceedings  paper  is  currently  only  available  from  the  conference 
website  and  its  availability  may  not  be  ensured,  the  paper  is  included  in  this  final  report.  Finally, 
recent  benchmark  model  comparisons  of  the  fast  ray  model  and  a  wave  propagation  model, 
based  on  the  trapped  modes  in  propagation  in  a  homogeneous  Pekeris  waveguide,  revealed  an 
significant  deficiency  in  the  implementation  of  fast  ray  model.  The  details  of  the  deficiency  and 
the  correction  to  the  implementation  of  the  fast  ray  model  are  given  in  a  section  that  follows  the 
proceedings  paper. 
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ENCLOSURE  2 


FAST  ACOUSTIC  RAY  MODEL  FOR  THE  SCATTERING  FROM  TARGETS 
IN  SYNTHETIC  APERTURE  SONAR  APPLICATIONS 
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Abstract:  Monostatic  synthetic  aperture  sonar  (SAS)  data  sets  for  the  scattering  from  a  target  in 
a  waveguide  can  be  simulated  via  a  fast  model  that  combines  an  acoustic  ray  approximation  for 
propagation  in  the  waveguide  with  far-field  scattering  from  a  target  in  free-space.  With  the 
assumption  that  acoustic  rays  arise  from  the  source  and  receiver  and  their  images,  the  problem 
of  wave  propagation  and  scattering  from  a  target  within  the  waveguide  is  replaced  by  a 
superposition  of  a  set  of free-field  scattering  problems.  Under  normal  operating  conditions,  the 
separation  distance  between  a  SAS  platform  and  a  target  is  large  compared  to  the  carrier 
wavelength  of  the  transmitted  signal,  and  hence,  the  free-field  scattered  signal  can  be  reduced  to 
a  far-field  approximation  where  a  spherically  diverging  wave  is  weighted  by  a  scattering 
amplitude.  The  scattering  amplitude  contains  all  of  the  information  about  the  scatterer  (e.g.,  its 
material  properties)  and  the  directionality  of  the  scattered  field.  A  scattering  amplitude  can  be 
obtained  from  a  direct  measurement  or  a  finite-element  analysis  of  a  target.  The  fast  ray  model 
allows  one,  who  is  interested  in  SAS  processing  algorithms,  to  generate  data  sets  with  variations 
in  the  environment  or  source-receiver-target  geometry  without  incurring  the  expense  associated 
with  the  collection  of  actual  SAS  data.  Results  from  the  fast  ray  model  will  be  compared  to  data 
and  finite-element  results  for  several  targets  in  shallow  water.  The  horizontal  range  to  the 
targets  and  time  gating  permit  the  air-water  interface  to  be  ignored. 


Keywords:  Acoustic  Ray  Model,  Synthetic  Aperture  Sonar,  Scattering  amplitude. 


1.  INTRODUCTION 


The  use  of  synthetic  aperture  sonar  (SAS)  has  become  a  standard  modality  for  the  detection  of 
targets  in  marine  environments  [1-3].  In  particular,  SAS  platforms  that  operate  at  frequencies 
below  50  kHz  and  with  a  wide  bandwidth  have  several  advantages  over  other  sonar  technologies 
(such  as  side-scan  sonar,  which  can  be  limited  in  range  and  range  resolution).  Low-frequency, 
broad-bandwidth  SAS  platforms  can  provide  longer  detection  ranges,  the  possibility  of  coupling 
into  elastic  modes  of  vibration  of  a  target,  and  deeper  penetration  into  sediments.  Longer 
detection  ranges  are  possible  due  to  the  lower  attenuation  of  low-frequency  sound  in  water.  This 
then  permits  rapid  wide-area  surveys.  In  addition,  the  range  resolution  of  a  SAS  platfonn  is 
related  to  the  bandwidth  of  the  transmitted  signal  and  is  independent  of  the  range  to  the  target. 
Excitation  of  the  elastic  modes  of  vibration  of  a  target  offers  a  mechanism  that  may  be  exploited 
by  classification  algorithms  to  differentiate  a  target  from  clutter.  With  greater  penetration  depths 
into  sediments,  the  detection  of  partially  and  completely  buried  targets  becomes  feasible.  Kargl 
et  al  [2]  recently  presented  a  model  to  investigate  the  response  of  buried  targets  to  a  sonar  signal. 
That  model  combined  first-order  perturbation  theory  for  scattering  from  a  rough  surface  and  an 
approximation  based  on  free-field  scattering.  This  paper  reports  on  a  model  that  combines  an 
acoustic  ray  model  and  free-field  scattering  for  targets  within  the  water  column. 

In  the  development  and  validation  of  automatic  target  recognition  (ATR)  algorithms,  access  to 
a  large  collection  of  SAS  data  sets  for  a  given  target  in  various  environments  and  various  targets 
in  the  same  environment  is  desirable  [4].  However,  the  in-situ  collection  of  data  in  a  marine 
environment  can  be  time  consuming,  cost  prohibitive,  infeasible  due  to  limited  access  to  a 
particular  environment,  and  subject  to  either  changes  to  the  environment  during  the  collection 
period  or  a  lack  of  measurement  of  environmental  conditions.  The  model  presented  below 
allows  one  to  fabricate  SAS  data  for  a  target  in  a  marine  environment  with  the  caveat  that 
passive  noise  and  reverberation  must  be  obtained  from  an  additional  model  (such  as  the  model 
described  in  [2]).  The  generated  SAS  data  sets  can  supplement  other  SAS  data  sets  and  aid  in 
the  development  and  validation  of  ATR  algorithms. 


2.  ACOUSTIC  RAY  MODEL  FOR  TARGET  SCATTERING 

When  the  wavelength  of  sound  is  much  smaller  than  the  depth  of  a  waveguide,  the  scattering 
of  sound  from  a  target  within  a  homogeneous  waveguide  can  be  approximated  by  an  acoustic  ray 
model.  The  waveguide  of  interest  here  is  a  homogeneous  layer  of  water  between  an  upper  semi¬ 
infinite  half-space  of  air  and  a  lower  semi-infinite  half-space  of  a  homogenous  sediment.  Figure 
1(a)  depicts  the  scattering  problem,  and  displays  the  direct  path  arrival  as  well  as  the  ray  paths 
that  interact  once  with  the  upper  and/or  lower  boundaries.  By  considering  image  sources  and 
receivers  reflected  about  the  boundaries,  one  can  associate  acoustic  ray  paths  with  these  images. 
Figure  1(b)  depicts  the  image  sources  and  receivers  and  their  associated  acoustic  rays.  Figure  1 
also  demonstrates  the  reduction  of  the  waveguide  scattering  problem  to  an  equivalent 
superposition  of  (possibly  infinitely  many)  free-field  scattering  problems.  In  the  model 
discussed  here,  the  homogeneous  sediment  has  been  modeled  as  either  an  attenuating  fluid  with  a 
frequency-independent  loss  parameter  or  as  a  fluid  described  by  an  effective  density  fluid  model 
[5]. 
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Figure  1:  (a)  Acoustic  ray  diagram  showing  the  direct  path  arrival  and  the  ray  paths  associated 
with  a  single  interaction  with  the  upper  and/or  lower  boundaries  of  the  waveguide,  (b)  Image 
sources  (Si,  S2,  ...)  and  receivers  (R ,,  If,  ...)  reduce  the  scattering  from  a  target  in  a  waveguide  into 
a  series  of  free-field  scattering  problems. 


For  an  image  source,  an  image  receiver,  and  a  target  located  at  ry,  and  rh  the  horizontal 
distances  are  Rti  =  \Rt  -  R,  |  and  Rj,  =  \R,  -  R,  |,  and  the  total  separation  distances  are  dti  =  \rt  -  r,  | 
and  dj,  =  \r,  -  r,  |,  respectively.  In  the  enumeration  of  image  sources  and  receivers,  i  =  0  and  j  =  0 
correspond  to  the  actual  source  and  receiver.  The  contribution  of  the  zth  source  and  the /h 
receiver  to  the  frequency  spectrum  of  the  total  scattered  signal  can  be  written  as 


P„(.a>)  = 


„i  exp(»at,,) 


Unu,Lm 


d 


jt 


jn(i)  rm(i)  CXp(/Vnt/;  ) 


unwr 


df, 


f{6ljn(plj,(o)Psrc{(o), 


(1) 


where  co  is  the  angular  frequency  and  Ps,x(o))  is  the  frequency  spectrum  of  the  transmitted  wave 
packet.  The  reflection  coefficients  at  the  upper  and  lower  boundaries  are  U{ 9)  and  L( 9).  Here,  0 
denotes  a  local  grazing  angle  at  the  boundary,  and  is  given  by  cos(0,)  =  Rtl  /  dti  or  cos(0y)  =  Rj,  / 
dj,.  The  m  and  n  exponents  indicate  the  number  of  interactions  a  ray  has  with  a  given  interface. 
The  time  delay  for  propagation  from  the  source  to  the  target  is  tn  =  dti  /  c  and  the  time  delay  from 
the  target  to  the  receiver  is  t,,  =  djt  /  c.  The  speed  of  sound  in  the  water  is  c.  In  (1),  the  sources 
and  receivers  are  assumed  to  be  point-like  with  omni-directional  directivity  patterns.  The 
scattering  process  is  a  convolution  of  a  free-field  scattering  amplitude /(0y,cj)y,co)  with  the 
incident  pressure  spectrum  Psrc( co).  Here,  0,y  and  <j);/  are  target-centered  angles,  which  are  related 
to  the  0,  and/or  0/  and  may  also  depend  on  the  orientation  of  the  target  within  the  waveguide. 

The  free-field  scattered  pressure  in  the  far- field  of  a  target  is  psrc  ~  pof(Qf)  exp (ikr)  /  r,  where 
the  wavenumber  is  k  =  co  /  c,/(0,cj))  is  a  scattering  amplitude  with  spherical  coordinate  polar  and 
azimuthal  angles  0  and  cj)  [2],  and  po  is  a  constant  carrying  the  units  of  pressure.  For  a  finite 
target,  the  scattering  amplitude  is  often  related  to  the  scattering  form  function  by/(0,(j))  =  aF 
(0,cj))  /  2  with  a  being  a  characteristic  dimension  of  the  target.  Inspection  of  psrc  shows  that  the 
target  is  point-like,  where  the  directionality  of  the  scattered  field  is  contained  within/(0,(j))  along 
with  infonnation  about  the  material  properties  of  the  target.  Finally, /(©^(ji^co)  can  be  obtained 
from  an  analytic  solution  to  a  scattering  problem  (e.g.,  scattering  from  a  spherical  target),  direct 
measurement  from  an  actual  target,  or  numerical  simulation  (e.g.,  finite-element  analysis). 


5 


3.  APPLICATION  OF  THE  MODEL  TO  CIRCULAR  SAS 


As  an  example,  the  waveguide  model  in  (1)  is  applied  to  a  monostatic  circular  SAS  (CSAS) 
scenario  where  the  duration  of  the  generated  time  signals  is  such  that  arrivals  from  the  air-water 
interface  are  excluded.  This  then  yields  four  rays  paths  where  f  =  0,1  and  j  =  0,1.  For  simplicity, 
we  use  do  =  dt o  =  dot  and  d\  =  dt\  =  d\t.  The  total  frequency  spectrum  is  then  given  by 


P{co)  = 


r  exp (ioatx)  ,  r,T(Q  w  exp (icatfi)  t2(Q  w  exp(iG*3) 
J\ - -p_ - +  ZL\Ug  )j2 — — —  +  L(dg)h 


dnd 


owi 


d{ 


Psrc(a), 


(2) 


with  t\  =  2do/c,  t2  =  (do  +  d\)/c,  and  t2  =  2d\/c,  and  the  grazing  angle  is  0„.  The  scattering 
amplitudes,  fk  =  //, ( 0 , <{) // , oa ) ,  depend  on  the  locations  of  the  sources,  receivers,  and  target,  and 
target  orientation  (in  the  Fraunhoffer  region,  and  equivalent  expression  with  some  useful  angle 
definitions  can  be  found  in  [6].  Eq.  3).  An  inverse  Fouier  transfonn  of  P( co)  then  gives  a 
generated  sonar  signal  that  includes  the  four  primary  acoustics  paths  for  a  target  near  an 
interface.  The  first  term  in  the  bracket  of  (2)  is  the  direct  path.  The  second  term  includes  the  two 
paths  that  interact  with  the  bottom  once.  These  paths  are  reciprocal  and  are  associated  with  a 
bistatic  scattering  direction.  The  last  term  is  a  backscattering  path  with  two  bottom  interactions. 

For  the  results  presented  here,  the  target  is  placed  at  the  center  of  the  CSAS  path.  It  is  a 
machined  aluminum  replica  of  a  solid  steel  artillery  shell,  which  was  deployed  during  Pond 
Experiment  2010  (PondExlO).  In  (2),  we  use  a  free-field  scattering  amplitude  derived  from  a 
hybrid  2D/3D  model.  The  hybrid  2D/3D  model  implemented  a  2D  finite-element  model  (FEM) 
prediction  of  the  target  response  on  a  surface  near  the  target,  and  then  propagates  the  target 
response  to  some  specified  range  using  a  discrete  fonn  of  the  3D  Helmholtz  integral  with  a  free- 
field  Green  function  [7-9].  The  hybrid  2D/3D  model  can  be  applied  to  targets  with  a  cylindrical 
symmetry.  A  target  drawing  and  the  FEM  mesh  is  depicted  in  Fig.  2.  This  scattered  pressure  is 
converted  to  a  scattering  amplitude 

Figure  2:  Engineering  rendering  of  the  aluminium 
replica  of  a  solid  steel  artillery  shell  and  the  FEM 
mesh.  The  nodes  of  the  mesh  are  concentrated 
near  the  fine-scale  structures  of  the  target.  The 
diameter  is  4  inches  and  the  length  is  16. 


and  placed  in  a  table,  f^,  with  0,  =  z'AO,  %  =_/A<|),  and  co*  =  kAco.  Given  the  source  and  receiver 
locations  along  the  circular  path  and  the  local  grazing  for  the  target  location,  values  of  0  and  (j) 
are  determined.  Linear  interpolation  was  then  employed  to  extract  scattering  amplitude  values 
from  the  table. 

Two  simulations  were  perfonned  with  the  radii  of  the  CSAS  path  of  5  and  10  m,  and  adjacent 
signals  were  recorded  at  A0  =  1°  along  the  circular  path.  The  source  and  receiver  were  co¬ 
located  3.8  m  above  the  sediment  and  were  assumed  to  be  omni-directional.  The  sound  speed 
and  density  of  the  water  are  c  =  1500  m/s  and  p  =  1000  kg/m  .  The  sediment  was  an  attenuating 
fluid  with  cs  =  1775  m/s,  pv  =  2000  kg/m  ,  and  8  =  0.008  for  the  loss  parameter.  The  source 
emitted  a  6  ms  chirp  with  a  carrier  frequency  of  16  kHz  and  bandwidth  of  30  kHz.  These 
parameters  roughly  correspond  to  experimental  conditions  during  PondExlO.  Figure  3  shows  the 
acoustics  color  templates  obtained  from  PondExlO  data,  FEM  simulations,  and  the  ray  model. 
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Acoustic  color  templates  depict  the  target  strength  as  a  function  of  frequency  and  a  target- 
centered  aspect  angle  projected  onto  the  water-sediment  interface  [6]. 

Inspection  of  Fig.  3  shows  that  both  the  FEM  and  ray  model  results  capture  much  of  the 
structure  observed  in  the  PondExlO  data.  This  structure  includes  both  geometrically  reflected 
acoustic  energy  and  acoustic  energy  that  couples  into  the  elastic  response  of  the  target. 
Comparison  of  Figure  3(b)  and  (c),  particularly  in  the  -10  to  -17  dB  range,  reveals  that  the  ray 
model  result  has  a  slight  smoothing.  This  is  a  consequence  of  the  discrete  representation  of/jjk 
and  the  linear  interpolation.  The  good  agreement  may  be  anticipated  because  a  free-field  FEM 
and  Helmholtz  integral  result  for  a  target  at  a  10  m  range  was  used  in  the  tabulation  of  the  /jjk. 


f  (kHz)  f  (kHz)  f  (kHz) 


Fig  3:  Acoustic  color  templates  for  the  aluminum  replica  at  a  10  m  range,  (a)  PondExlO 
measurement,  (b)  FEM  and  Helmholtz  integral  result,  (c)  Ray  model  using  the  tabulated 
form  function.  The  tail  of  the  target  is  pointing  at  the  source/receiver  location  at  90°;  while 
the  nose  points  at  the  source/receiver  at  270°. 
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Fig.  4:  Acoustic  color  templates  for  an  aluminum  replica  at  a  5  m  range,  (a)  PondExlO 
data,  (b)  FEM  model  and  Helmholtz  integral  result,  (c)  Ray  model  using  the  tabulated 
form  function. 


Figure  4  shows  the  PondExlO  data  for  the  aluminum  replica,  FEM  and  Kirchoff-Helmholtz 
integral  simulation,  and  the  ray  model  given  in  (2)  for  the  target  at  a  5  m  range.  An  important 
aspect  of  the  ray  model  simulation  is  that  it  used  the  tableyijk  that  was  derived  from  a  free-field 
FEM  and  Helmholtz  integral  simulation  for  a  target  at  a  10  m  range.  This  is  possible  because  the 
scattering  amplitude  is  range  independent  and  contains  only  infonnation  about  the  directionality 
of  the  scattering  and  material  properties  of  the  target.  The  construction  of  the  model  in  Eq.  (2) 
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explicitly  accounts  for  interactions  with  the  water-sediment  interface.  Thus,  only  a  single  free- 
held  measurement  from  a  target  or  a  single  free-held  FEM  and  Helmholtz  integral  simulation  is 
sufficient  to  provide  a  scattering  amplitude  that  can  be  used  to  rapidly  predict  the  scattering 
from  a  target  near  a  bound  at  any  range  (provided  the  range  meets  the  far  field  requirement). 


4.  CONCLUSIONS 

Section  2  presented  a  model  that  combines  a  ray-based  description  of  wave  propagation  in  a 
waveguide  with  a  far-field  approximation  for  the  free-field  scattering  from  a  target.  In  Sec.  3, 
the  model  given  in  Sec.  2  was  applied  to  a  circular  SAS  example  for  a  proud  target  where 
interactions  with  the  upper  boundary  could  be  ignored  (e.g.,  due  to  time  gating).  Figures  3  and  4 
demonstrate  that  the  ray  model  captures  the  structure  observed  in  both  experimental  data  and  an 
alternate  simulation  that  used  the  hybrid  2D/3D  model.  At  the  time  that  Fig.  3(b)  and  4(b)  were 
computed,  a  hybrid  2D/3D  model  result  took  2  to  3  days  to  complete.  The  ray  model  results 
depicted  in  Fig.  3(c)  and  4(c)  took  on  the  order  of  30  CPU  seconds.  Thus,  the  ray  model  retains 
most  of  the  high  fidelity  of  the  hybrid  2D/3D  model  at  a  fraction  of  the  computational  resources. 
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Enclosure  3 


Fast  Ray  Model  Deficiency: 


A  fast  ray  model  developed  in  FY 1 1  assumed  a  homogeneous  layer  of  water.  Rays  are 
assumed  to  travel  in  straight  line  segments  from  a  source  to  a  receiver.  For  an  interaction  of  a 
ray  with  an  interface,  an  appropriate  reflection  coefficient  is  included  in  the  received  signal.  For 
an  image  located  at  r,  and  a  receiver  located  at  #y,  the  horizontal  and  total  separation  distances 
are  | R,  -  R,\  and  |#y  -  r,|,  respectively.  The  contribution  of  the  /th  image  to  the  frequency 
spectrum  at  the  receiver  can  be  written  as 


Pi(oj)  =  Psrc(aj)s0 


An(Gi)Bm(0i ) 


exp  (iaX,) 


(1) 


where  Psrc( co)  is  the  frequency  spectrum  of  the  transmitted  wave  packet  and  so  is  a  reference 
distance  used  in  the  calibration  of  the  source  level.  The  reflection  coefficients  at  the  upper  and 
lower  boundaries  are  A(Q,)  and  f?(0,),  where  the  local  grazing  angle  8,  at  a  boundary  is  given  by 
cos(0,)  =  | R,  -  Rj\  /  |#y  -  r,\,  and  the  m  and  n  exponents  indicate  the  number  of  interactions  a  ray 
has  with  a  boundary.  The  time  delay  for  propagation  is  U  =  \rr  -  r,\/c  with  c  being  the  speed  of 
sound  in  the  water. 

During  FY13,  simulations  of  scattering  from  a  15-m  radius  fluid  sphere  in  a  4000-m  deep 
Pekeris  waveguide  with  source-receiver-target  separation  distances  on  the  order  of  100  km 
revealed  the  implementation  of  the  propagation  component  of  the  fast  ray  model  violated 
reciprocity  under  certain  conditions.  If  the  source  and  receiver  were  located  below  the  mid¬ 
depth  of  the  Pekeris  waveguide,  the  implementation  satisfied  reciprocity.  If,  however,  either  the 
source  or  receiver  or  both  were  above  the  mid-depth  of  the  waveguide,  then  reciprocity  failed. 

The  failure  was  caused  by  a  change  in  the  enumeration  of  image  sources,  the  method  used  to 
detennine  the  reflection  coefficient  exponents,  and  the  selection  of  images  that  contribute  within 
a  time  window  of  interest.  The  depth  of  the  waveguide  is  h,  and  the  actual  source  is  located  at  zs 
(>  0)  with  the  seafloor  at  z  =  0.  The  first  image  was  below  the  seafloor  at  z\  =  -  zs.  Additional 
images  were  added  in  quartets,  because  images  planes  were  added  in  pairs,  symmetrically  about 
the  seafloor.  The  gth  quartet  of  image  was  added  to  the  set  of  image  with  the  following  z- 
coordinates:  2 qh  +  z„  2 qh  -  zs,  -2 qh  -  zs,  and  ~2qh  +  zs  with  q  =  1 ,  2,  3,  ....  The  ordering 
assumed  that  zs  <  h/2,  and  the  detennination  of  the  reflection  coefficients  assumed  that  the  ray 
associated  with  the  first  image  always  interacted  with  the  seafloor  and  the  rays  in  the 
enumeration  of  the  remaining  image  source  then  alternated  with  interacting  with  the  air-water 
interface  and  seafloor.  As  long  as  zs  <  hi 2,  z\=  —  zs,  and  the  addition  of  quartets  followed  the 
prescribed  ordering,  then  the  images  were  sorted  in  a  manner  consistent  with  the  detennination 
of  the  m  and  n  exponents  in  Eq.  (1),  and  reciprocity  was  ensured. 

In  the  FY13  simulations  with  separation  distances  on  the  order  of  100  km,  it  was 
advantageous  to  consider  only  a  finite  length  time  window,  tw,  with  a  defined  time  offset,  t0.  It 
then  followed  that  only  images  with  propagation  delays  f  =  | #y  —  r,|/c  within  the  time  window 
needed  to  be  considered;  and  so,  the  set  of  image  were  sorted  such  that  tQ  <  t,  <  tw  and  t,  <  ti+\ 
When  zs  >  h/2,  the  new  sorting  led  to  a  first  image  with  zi  =  2 h  -  zs,  and  the  ray  associated  with  it 
reflects  initially  from  the  air-water  interface.  The  implementation  of  the  new  model  failed  to 
updated  the  determination  of  m  and  n,  and  thus,  reciprocity  was  violated. 
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The  set  of  images,  which  contribute  to  the  pressure  spectrum,  are  selected  from  the  set 
constructed  in  the  following  manner.  First,  for  a  source  located  at  rs  =  (xs,ys,zs)  and  0  <zs<  h, 
coordinates  for  all  images  have  xt  =  xs  and  v;-  =  ys.  The  coordinates  for  N  images  (including  the 
actual  source)  are  first  constructed  in  an  unsorted  manner: 

Zj  =  zs ,  (actual  source) 


V3  =2qh  +  zs, 
zq+4  =  2qh  —  z  s , 

V  s=<2qh-zs), 
V6  =~(2<lh  +  zs)’ 


(3) 


where  q  =  0,  4,  8,  ...  and  q  +  6  <  N.  Next,  the  images  are  sorted  by  the  d,  =  \z,  -  z i|  with  i  =  2,  3, 
...,  N  and  di  <  dm.  This  then  guarantees  the  ordering  of  delay  times:  ti  <  ti+ 1. 

The  detemination  of  the  reflection  coefficient  exponents  then  follows  from  the  symmetry 
imposed  by  the  seafloor  at  z  =  0  and  whether  Z2  <  zs  or  Z2  >  zs.  If  the  former,  the  ray  tracing  of 
the  first  few  rays  yields  the  sequences  of  reflection  coefficients: 


i  =  1, 

(actual  source) 

i  =  2. 

B, 

i  =  3, 

A, 

i  =  4, 

AB, 

i  =  5, 

BA, 

i  =  6, 

BAB , 

i=7. 

ABA, 

i  =  8, 

ABAB, 

i  =  9, 

BABA, 

i  =  10, 

BABAB , 

i  =  ll. 

ABABA, 

the  latter  situation,  the  ray  tracing  gives 

i  =  1, 

(actual  source) 

i  =  3, 

B , 

i  =  2. 

A, 

i  =  5, 

AB, 

i  =  4, 

BA, 

i  =  l. 

BAB, 

=  6, 

ABA, 

i  =  9, 

ABAB, 

i  =  8, 

BABA, 

/  =  11, 

BABAB, 

i  =  10, 

ABABA, 

(4) 


The  violation  of  reciprocity  in  the  benchmarks  was  caused  by  an  unconditional  use  of  (4)  where 
the  location  of  the  actual  source  within  the  waveguide  was  ignored. 
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